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Abstract 

^ . We present classes of discrete reversible systems which are at the same time chaotic 

I and solvable. 

0^ ' 
O 

o 

^! 1 Introduction 

Chaos and solvability are antithetical notions and their coexistence in a dy- 
namical system may sound paradoxical. The term chaos is traditionally used 
I to designate systems which exhibit crucial dependence on the initial conditions 

j_j ■ manifested through exponentially diverging trajectories and extreme instabil- 



ity [1]. On the other hand, solvability is often associated to smooth, regular 
behaviour, related to the existence of invariants, and is usually coming from 
integrability. 

In what follows we shall show that the explicit solvability of a mapping is 
not incompatible with a chaotic behaviour. This is done in the same spirit 
as in [2,3], but we will give example of reversible systems (i.e. there exists 
a similarity transformation between the forward evolution and the backward 
evolution), exhibiting at the same time features of chaos (e.g. sensitive de- 
pendence on the initial conditions, positive algebraic entropy, ergodicity), and 
which are solvable. Our point is neither to recall the possible coexistence of 
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KAM tori and chaotic regions which is well known [4] , nor to give the ultimate 
definition of integrability, solvability or chaos. 

Our point is to present reversible systems which lie on the border of solvabil- 
ity/integr ability and chaos. 



2 Specific examples 

We will use maps which arc constructed from recurrences, that is to say se- 
quences where each term is given as a function of the previous ones: 

If Xn+i is a function of only a;„, we have a one dimensional map, also called a 
system of order one. li u — 1 (resp. v — r) we may define from (1) a map in 
two dimensions (resp. r + 1 dimensions) by 

\Xini X11—I1 ■ ■ ■ 1 Xn—u\ ^ [^^n+l) ■^ni ■ ■ ■ i ■^n—v+l\ (2) 

The space of initial conditions is of dimension u + 1. 
2.1 One dimensional maps 



We start with a well-known one dimensional example to illustrate the fact 
that solvability is compatible with chaotic behaviour: 



Xn+\ — 2 — 1 (3) 

This map is know to be chaotic [5]. Using the similarity between (3) and the 
doubling rule for cosine: cos( 2 ui) =2 cos^ a; — 1, the solution of (3) is given by 
Xn — cos2"q; where a is some constant determined by the initial conditions. 
We shall not dwell upon the exponentially fast loss of the memory of the initial 
conditions due to the presence of the 2" factor: the existing literature covers 
the topic in an exhaustive way [1]. It remains that we have here an example 
of a system that is explicitly solvable and which has a chaotic behaviour. 

The previous example is not an isolated occurrence. Whole families of first- 
order (one-dimensional) systems exist which are solvable while at the same 
time exhibiting exponentially fast loss of the memory of the initial conditions. 
There exist a number of results in this direction for maps of the interval [6,2,3]. 
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It was shown that the only nonhnear polynomial maps where the solution can 
be explicitly given are the map Xn+i — x"^ and (3), up to a homographic 
transformation. 

For rational maps his classification is based on results of Ritt [7] . For degree 
one, the only solution is the homographic map. For degree 2, there exist 8 
different recurrences, the simplest ones being 



Xn+l^^{Xn-—) (4) 
Zl Xn 

1 1 

Xn+l^-j{Xn-'2^ ) (5) 

4 Xfi 

Their solutions can be given in terms of the WeierstraB elliptic function p(^), 
defined through the equation [ip'Y = Ap^ — Ap. Equation (4) has solution 
— p((l + i)"z). Similarly the solution of (5) is Xn — + i)"'z), and so 

on. 



The major problem with the above mappings is that they are not invertible. 
While one can define the image of a given point in an unambiguous way, the 
same is not true for the inverse evolution. Moreover the number of preimages of 
a point grows exponentially fast (a property which was deemed incompatible 
with integrability [8]). In the next paragraph we will exhibit reversible maps, 
defined by birational transformations showing features of chaos and which are 
solvable. One condition to construct such maps is to consider transformations 
in more than one dimension. 



2.2 Two-dimensional maps 



In the case of transformation (3), the key ingredient was the doubling relation 
LUn+i = 2a;ji (or a higher multiple uJn+i = kuJn)- We start from a linear equation: 



LUn+l + ^n-\ = kuJn (6) 

The solution of (6) is straightforward: a;„ = aA" + 6A" where A± = (/c ± 
VA;^ — 4)/2. Exponentiating (6) and setting x — we obtain the recurrence: 

•^n+l-^n— 1 •^n C^) 



The chaotic character of this map (for k > 2) can be assessed easily through 
the computation of its algebraic entropy [9]. This quantity is a measure of the 
complexity of the map and is given by e = lim„^oo log ((i„)/n where d„ is the 
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degree of the n-th iterate. In the case at hand we find that e = log(pi), with 
the pi the larger of the two roots defined above. This leads to e > for k > 2. 
In fact the same value for the entropy is obtained for all the maps derived 
from (6). 

While map (7) is rather trivial, it is possible to construct a much more inter- 
esting one by setting 

Xn = tanun (8) 



We get the recurrence: 

= fkiXn) (9) 



where fk is a rational function of Xn depending on the value of A;. It is simply 
the expression of tanA;a;„ in terms of tanti;„ = Xn- For the first few values of 
k we have fi{xn) = Xn, f2{xn) = 2xn/(l - xl), h{xn) = (3xn - x^) / {I - 3x^) 
etc. The case k — 1 corresponds to a trivial map which is periodic with 
period 3. The case k = 2 can be easily integrated. There exists an invariant 
c = (xn-i—Xn) / (l+Xn-iXn) , which allows to reducc the map to a homographic 
one and solve it completely. 

The case A; = 3 is more interesting. Indeed the map constructed from the 
recurrence 



Xn+l 



SXfi x^ 3Xjj) 
1 3a;^ + {^Xn xf^) Xn—l 



(10) 



is both chaotic and solvable (since it is linearizable). The generating function 
of the sequence of degrees [10] may easily be inferred from the first terms of 
this sequence: 



^~ (1-S)(1 + S + S2)(1-3S + S2) ^ ' 

The value of the algebraic entropy e = log((3 + -\/5)/2) is read off from (11). 
This positive entropy is a sign of chaotic behaviour. 

Similar maps can be constructed for /c > 3. They will have a positive algebraic 
entropy and be linearizable. 
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3 Graphical analysis 



Figure 1 shows the orbit of an arbitrary initial point under 10^ iterations of 
the evolution (10), plotting the pairs {xn,Xn+i) in the two-dimensional plane. 
It shows that the orbits tend to fill phase space. 



3.0 




3 

Figure 1: A typical orbit of a point under map (10) 



Figure 2 shows the images of a segment under a dozen of iterations of (10). 
Pushing the iteration further does not change the qualitative features of the 
image. It just increases the density of lines. 
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The fact that the map (10) has an algebraic entropy of log((3 + -\/5)/2) invites 
us to compare it to the map described in [11], which has the same algebraic 
entropy, and is chaotic: 



X, 



n-1 



a 



[12) 



with a a constant. 
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The typical orbit of an arbitrary point was shown in [11], presenting large 
chaotic regions, and we do not reproduce it here. 

Figure 3 shows the images of a segment of a straight line under 27 iterations 
of (12), and is to be compared with Figure 2. The overall pictures differ quali- 
tatively. A detailed analysis of the behaviour of the iterates shows that in the 
case of the map (12) the evolution is slow till the points get near one of the 
singularities, in which case they get a major boost, which results in the rich 
structure of Figure 3. On the contrary, in the case of map (10) neighboring 
points are uniformly repelled and the overall results are regular and smooth. 
Both have a strong dependence on the initial conditions, but in different ways. 

All numerical calculations leading to Figures 2 and 3 have been performed 
with multi-precision [12] arithmetics so as to guarantee the reliability of the 
results, and a test was performed on all image points: after an evolution of N 
steps, we "go back in time" the same number of steps to recover the initial 
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data. We adjusted the precision in such a way that the true initial data and 
the result of this round trip did not differ by more than 10~^. 

The distinguishing feature of (10) compared to (12) if that (8) has an integral 

(a;„+i - A+ a;„)(a;„+i - A_ a;„) = est (13) 

with A± = (3 ± v^)/2. The c<j-plane is foliated by invariant curves, but the 
picture is scrambled by the transformation (8) since the tangent function is 
periodic. It is easy to understand the aspect of Figure 1. The curve (13) is 
transformed by (8) into 

(arctanx„_^i — A+ arctan2;„)(arctanx„+i — A_ arctanx„) = k (14) 

with Xn+i and Xn the two coordinates of the plane of Figure 1, and k a constant. 
Such a curve cuts a line a;„ = ^ at an infinite dense set of points. 

One may notice that there exist two holomorphic foliations which are left 
invariant by (9), as in the analysis of [13], and our construction exemplifies 
their result. 

One interesting feature of transformation (8) is that it is not a mere change 
of coordinates. It takes a birational map into a birational one, changing the 
algebraic entropy. 



4 Arithmetical analysis 

The difference between maps (10) and (12) can be illustrated by an analysis 
based on the approach recently introduced by Roberts and Vivaldi [14]. These 
authors have studied the effect of the existence of rational integrals of motion 
for rational maps when the evolution is considered over a finite field. The 
simplest realization of such an evolution is through integer arithmetics modulo 
some prime integer p. The basic observation of [14] (see also [15]) is that if 
there exists a rational invariant, the orbit is confined to an algebraic curve, 
and the genus g of this curve is at most g < 1, if the original map is of infinite 
order. Such curves over finite fields have a maximum number of points (the 
Hasse-Weil bound HW{p,g) = p + 1 + 2 g y/p), and as a consequence the 
number of points on the same orbit is also bounded by HW{p,g). In short: 
algebraically integrable maps have a large number of orbits and those arc 
short. Chaotic maps have a smaller number of orbits and they are longer. 
Notice that this fits with the idea that chaotic orbits may explore the whole 
phase space, contrary to what happens in the integrable case. 
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We have performed a sampling of initial points for increasing values of p and 
plotted the mean value of the length of orbits (with the rule to terminate the 
iteration when meeting a singular point or closing a loop), for the maps (10), 
(12), and one additional map which is known to be algebraically integrable 
(so called McMillan map), given by the recurrence: 




Figure 4: Mean length of orbits vs p 



Figure 4 shows the values of the mean length for increasing p, normalized 
by dividing by HW{p, 1) = p + 1 + 2 y/p, for maps (10), (12), and (15). 
It discriminates between (10) and (12) by showing that the former verifies a 
virtual Hasse-Weil bound while the latter does not. The curve corresponding to 
map (15) is there for reference, as well as the line corresponding to a constant 
normalized length of 1. Clearly, although the orbits are not confined to any 
invariant algebraic curve, the arithmetic test of [14] places the map (10) in 
the integrable class together with (15), in contrast to (12). 
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A more detailed analysis of the statistics of the length of orbits will be done 
elsewhere. 



5 Generalizations 



We can construct generalizations of (9), based on the properties of elliptic 
functions. Introduce Xn = p(cc'„) and = p'iujn) where p is the WeierstraB 
elliptic function with elliptic invariants g2 and g^. The addition formulae for 
the WeierstraB eUiptic function yield: 



1 / hk{xn,yn) +y n—1 \ /. / \ 

4 \Jk{Xn,yn) - Xn-lJ 

Xn+l{hk{Xn,yn) +yn-l) 

Vn+l = —FT N (16) 

•^n— 1 JkKp^m yn) 

fk{Xn,yn)yn-l + 

fkipni yn) 

where the functions fk and hk are the expressions of p{ku!n) and p'{ku!n) in 
terms of Xn = p{oUn) and yn = p'{oUn)- 

We have to ensure that both {xn,yn) and Xn-i,yn-i he on the same elliptic 
curve i.e. 



yl = ^xl- g2Xn- 93 (17) 
yl-i = 4a;^_i - g2Xn-i - 93 (18) 

Equations (17,18) together with the formula for the duplication, tripUcation, ... 
of the arguments yield the value of the fk and hk- Iteration (16) defines a map 
in four variables which has two algebraic invariants given by solving (17,18) 
in terms of g2,g3' 



92 = + 4{xl_, + XnXn-l + 4) (19) 

2 2 

93 = z 4x„x„_i (X„_i + Xn) (20) 



For k = 2 we have /2 = —2x + z'^/{Ay^) and /i2 = —y + 3xz/y — z'^/{4y^) where 
z stands for p"{u!n) = 6p^(a;„) — 92/2 = Qx^ — g2/2. In this case there exists 
an additional invariant 
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C= f ^"'_^^^" )'-4(x„ + x„^i) (21) 

and the map is integrable with vanishing entropy (quadratic growth of the 
degree). 

For = 3 we have /s = x + 4y'^{12xy'^z — Ay'^ — z^)/{12xy'^ — z^)'^ and /13 = 
—y — A.yiVlxy^z — 8y^ — z^){12xy'^z — Ay'^ — z^)/{12xy'^ — z'^Y- This case, and 
actually all cases with /c > 3, have positive entropy and solvability is ensured 
through the relation to (6). 



6 Conclusion 



We have shown that there exist infinite families of rational maps which, at the 
same time, have positive algebraic entropy, present features of chaos, and are 
solvable. Their solvability is related to a reduction to a linear equation through 
the appropriate non rational transformations, but they remain reversible. 

While the examples we exhibited here are based on specific Ansatze, there 
exist infinite families of solvable mappings with positive algebraic entropy. As 
a matter of fact one could perform the same derivation using any function for 
which one can express f{x + y) in terms of f{x) and f{y)- 

Open questions remain, like what is the meaning of the statistics of the length 
of orbits in the arithmetic approach of [14,15]. We will return to that in some 
future pubhcation. 

Acknowledgments: We acknowledge stimulating discussions with C. Favre. 
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